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1 Introducntion 



Determinant representations of correlation functions of one dimensional quantum 
integrable models were studied extensively 0. Especially Essler et. al. studied the 
correlation functions of XXZ model with periodic boundary conditions M. In the 
calculation of scalar products of the XXZ model, the relation between the scalar 
products of the model and the partition function of the six- vertex with domain wall 
boundary conditions is essential 0. Izergin et. al. used this relation to obtain the 
determinant formula for the partition function of the six-vertex model with domain 
wall boundary conditions ||. 

On the other hand, integrable models with open boundary conditions are inter- 
esting. Integrable model with open boundary conditions are related to the Bjy type 
Weyl group || and in condenced matter physics they are related to the impurity 
problems in the Luttinger liquid 0. 

Correlation function of the 1 + 1 dimensional delta function interacting Bose 
gas with open boundary conditions were studied in ||. 

For the XXZ model with open boundary conditions, the algebraic Bethe ansatz of 
the model can be used to calculate the partition function of the six- vertex model with 
reflecting boundary condition at one boundary and domain wall boundary conditions 
at other three boundaries. (In this paper we call six vertex model with this type of 
boundary conditions six vertex model with reflecting end.) In this paper, we obtain 
the determinant formula for the above six-vertex model using above relation. 



Quantum Inverse Scattering Method for the Open 
XXZ Model 



Let us start with a brief review of the Quantum Inverse Scattering Method for the 
XXZ model with open boundary conditions ||. The trigonometrical solution of the 
Yang-Baxter equation is given by 
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where i?oo' (/•*) acts as a linear operator on the tensor product of two auxiliary spaces. 
The L operators are given by L n (X) = smh.(X)R 0n (X)P 0n where denotes the auxil- 
iary space, n denote the n-th quantum space and 
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Pon = (1/2) (lo ® In + o"o ® a n + a o ® °"n + °"o ® U n) ls a permutation operator. 
Thus 

L n {\) = -{[sinh(A-2i77)+sinh(A)]lo®l„-sinh(2i77)[cr^®< + (7^®(7^ 

+ [sinh(A - 2irj) - sinh(A)] o% ® <} . (3) 

As a consequence of the Yang-Baxter equation, the L operators satisfy the fun- 
damental commutation relations 

R{\-n)(L n {\)®LM) = (L n (ji)®L n (\))R(\-n). (4) 

The monodromy matrix of inhomogeneous XXZ model is constructed from the 
L operators as: 

T{\,u) = L L {X-u L )L L - 1 {X-u L - 1 )---L 1 (X-u 1 ). (5) 

This matrix satisfies fundamental commutation relations similar to (f|) 

R(X — n) (T(A, u)) ® T(n, u;)) = (T(ji,w)®T(\,w)R(\-n). (6) 

For the model with open boundary conditions, in addition to the Yang-Baxter 
equation, the reflection equation || plays an important role, 

R 12 (X- v)K 1 (X)R 12 (X + riMfi) = K 1 ( f i)R 12 (X + fi)K 1 (X)R 12 (X-fi), (7) 

where Kj(X) is the 2x2 matrix which act to the j-th auxiliary space and is related 
to the reflection of the particle with rapidity A at the boundaries. In the present 
paper, we concentrate to the diagonal solution of the reflection equation 



K(X,0 



sinh(£ + A) 



sinh(f) V sinh (£ - A ) / ' 

For the model with open boundary conditions, Sklyanin's monodromy matrix 



(9) 



where 



T(A, u) — L a (A + ooi)L 2 (X + cu 2 ) ■ ■ ■ L L (X + uo L ) 

xT-\-X,u) (10) 

is fundamental, where A(X, u>), B(X, u>),C(X, U)),T>(X, u>) are operators which act on 
the quantum spaces. 
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U(\,uj) also satisfies reflection equation 

Rn(X - /^)C/i(A, u)fl 12 (A + //)C/i(^ w) 

= ^(/i, u>)i2i 2 (A + //)[/! (A, u;) J R 12 (A - /i). (11) 

The reflection equation is rewritten as the commutation relations of the matrix 
elements A, B, C, V, 

[A(X),A(»)) = -[V(X),V(ri] 

smh(2ir]) 



sinh(A + fx — 2ir)) 
sinh(2ir]) 



{B{\)C(p.) - B(im)C(\)} 
{C(\)B(ri-C(n)B(\)} 



sinh(A + fx — 2irf) 

= sinh ( A - V) sinh ( A + V) \A{\) V(u)\ (12) 
sinh 2 (2zr]) — sinh(A + fi — 2%r\) ' 

[B(A),B(aO] = [C(A),C(ai)] = (13) 

sinh(2z?]) 

[i3(A) ' C(/i)J = sinh(A - S sinh(A + Ji - 2mj) X 
x{sinh(A + //)[£>(A).4(//) - V(fi)A(X)\ 

+ sinh(A - //)[^(A)^(//) - V(fi)V(X)]} (14) 

(15) 

From the above commutation relations (|T3D Sklyanin's transfer matrices 

r(A) = troCKX-A + 2irj)U(X)) (16) 
= (1/ sinh(£ + )) [sinh(£ + - A + 2ir])A(X) + sinh(£ + + A - 2ir])V(X)} (17) 

are commutative for any values of spectral parameter [r ( A), r (//)]. 

Consequently, r(A) generate the infinite number of conserved quantities. For 
example the Hamiltonian is given by 

^(A)| fc „ 

1 L ^ 

= - sinh(2 ^) + ^n+i + cosh(2^)« +1 } 

+ coth(£_)<7i - coth(£ + )(x2 + constant (18) 

In the algebraic Bethe ansatz, we construct the simultaneous eigenstates of the 
conserved quantities by applying the some elements (creation operator) of the mon- 
odromy matrix to the pseudovacuum. We choose as a pseudovacuum the state whose 
spins are all up 

|0) = ®t=il+)r ( 19 ) 
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This pseudovacuum is the eigenstate of A and V and is annihilated by C 



„4(A)|0) = a(A)|0>, 
2>(A)|0) = d(A)|0), 

C(A)|0> = 0, (20) 

where 

L 

a(A) = sinh(£_ + A) {sinh(A — uj n — 2ii]) sinh(A + uj n — 

n=l 
L 

d(X) = sinh(£_ - A) J] {sinh(A - u n ) sinh(A + uj n )}. (21) 

n=l 

In other words B has the role of the creation operator and the eigenstates of the 
Hamiltonian are of the form IT^Li £>(^n? <•<-') |0) if rapidities A n satisfy the Bethe ansatz 
equations. 

In the algebraic Bethe ansatz eigenstates and dual states of eigenstates are re- 
spectively given by 

B(X L ) ■ ■ -B(Ai)|0) (22) 

and 

(0\C( f i 1 )---C^ L ). (23) 



3 Inhomogeneous six-vertex model with reflect- 
ing end 

For the scalar products between the eigenstates the following relation holds 

(0\C(^)---C( f x L )B(X L )---B(X 1 )\0) 

= (0|C(/xi) • • .C(^)|0>(0|B(A L ) • • .B(A0|0> (24) 

where |0) is the quantum state with all spins are down Yin=i \ ~ )• Note that 

Z L ({\},{u>}) = (Q\B(\ L )---B(\ 1 )\0), (25) 

is nothing but the partition function of the inhomogeneous 6 vertex model on a 
L x 2L lattice with domain wall boundary condition for the lower, upper and left 
boundaries and reflecting end for the right boundary (see Fig. 1.). 
S(An) =° n (+M\u)\-)° n where |+)° and |-)° 
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are the n-th auxiliary states with spins + and respectively — (°(+| and °(— | are 
dual states with spins + and respectively — ). Here we assign the auxiliary space to 
each B(X n ). 

Now we derive the recursion relations to the partition function Zl({\}, {uj}). 
If the inhomogenuity ojl is choosen to be ojl = —\\ then the partition function 
of the Lx2L lattice is related to the partiton function of the L—lx 2(L — 1) lattice. 
The L operators are regular i.e. 

L L (0) = -smh(2ir))P oL . (26) 

Thus 

L 

#(Ai) n i+ >3 u=-\i 

3=1 

L 

= - sinh(2^)°(+|L L (A 1 - u L ) ■ ■■L L _ 1 (\ 1 + c^P^-)? J[ \+)j 

3=1 
L-l 

= - sinh(2^)°(+|L L (A 1 - u L ) ■ ■ -L L -i(Ai + |+> 

3=1 



3- 

(27) 



Noting that |+)i|+)„ is the eigenstate of L n (\i) 

^n(Ai)|+)i|+>„ = sinh(A! - 2i77)|+)i|+) n , (28) 

and 

i(+|L L (A n - u L )\-) L = 2sinh(A n - (30) 

We obtain 

L 

n i+ >3 u=-xi 

3=1 

sinh(^i + Ai) . . 

■ smh(2^) smh(2Ai) 



sinh(£_) 

L-l 

x TT sinh(Ai — ujj — 2ir\) sinh(Ai + ujj — 2irj) 

3=1 

L-l 



X 



-) L n i+>* (si) 

3=1 
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Using 

L L (X n + u L )\—) n \-) L = sinh(A n + u L - 2ir))\-) n \-) L , (32) 
together with equation fl3"0|) and the locality of the L operators we find the relation 

= sinh(A„ + u L - 2irj) sinh(A„ - uj L )\-) L B L -i(X n ), (33) 

where 

B L -i(X n ) = n (+\LL-l(\n-VL-l)'--L 1 (\ n -Ui)K(\ n ,£-) 

xL!(A n + u{) ■ ■ • L L _ x (A n + w x _i)|-) n . (34) 

lemma 1-a 

The partition function of the model on the L x 2L lattice is related to hte one 
on L — 1 x 2(L — 1) lattice 



^({A},M)U = _ Al 

sinh(£_ + X L 



. sinh(2w/) sinh(2A L ) 
smh(£_) 

i-i 

x Y[ sinh(Ai — Uj — 2ir)) sinh(Ai + ujj — 2ir\) 

3=1 

L-l 

x Y[ sinh(A n + u L - 2irf) sinh(A n - u L )Z L -i({X}, {u}). (35) 

n=l 

Here Zm({X}, {uj}) is the partition function for the model on M x 2M lattice 

M M M 

Z M ({A},M) = IL'H II Bu(*»M II (36) 

3=1 n=l i=l 

which depend on Ai ~ Xm and uj\ ~ cj A /. 
lemma 1-b 

Folowing the same line of argument as above we have other type of recursion 
relations 
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^l({A},M)U =A: 

sinh(£_ - X L ) 



sinh(2i^) sinh(2Ai) 



sinh(£_) 

L-l 

x Y[ sinh(Az, — uij — 2ir/) sinh(AL + tOj — 2ir\) 

3=1 
L-l 

Y[ sinh(A„ - u L - 2irj) sinh(A n + ^ L )Z L _i({A}, {a;}). (37) 

n=l 

Lemma 2 

Z 1 ({A},M) = 1 (0|S(A)|0) 1 (38) 

is explicitly calculated as 

Zi({A}, {u}) = - SmI " (2 ^ [smli(e + A) sinh(A - w) + sinh(£ - A) sinh(A + jffl) 
Lemma 3 

(0\B(Xl) ■ ■ •B(Ai)IO) depend on each A n as 

(0|S(A L )---S(A 1 )|0) 

= e~ 2LXn P 2L (e 2Xn ) (40) 

where P2l(x) is the polynomial of degree 2L. 

This is proved along the same line of argument as in reference ||. 

4 Determinant formula for the six vertex model 
with reflecting end 

The lemmas in the last section determine the partition function of the model uniquely. 
Theorem: The partition function of the six vertex model with reflecting end is 
given by 

Zl({\},{"}) 

L L 

J ! j J sinh(Aj — u n — 2irj) sinh(A 3 - — u n ) sinh(Aj + u n — 2ir\) sinh(Aj + u n ) 

3=1 n=l 

x[Y{ sinh(Aj - X k ) sinh(Aj + X k - 2ir\) J\ sinh(u;„ - uj m ) sinh(u; n + uj m )}~ 1 

j>k n>m 

xdetM (41) 
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where 



M = sinh(2^) 

Jn sinh(£) sinh(Aj — u n — 2irj) sinh(Aj + u n — 2irj) 
{ sinh(£_ + Xj) sinh(£_ — Xj) 
\ sinh(Aj + uj n ) sinh(Aj — uj n ) 

5 Conclusion 

In this paper the determinant formula for the partition function of the six-vertex 
model with reflecting end is given. It is the determinant formula corresponding to 
the Quantum Inverse Scattering Method with open boundary conditions. 

In the model with open boundary conditions, as opposed to the model with 
periodic boundary conditions, one can not determine the scalar product from the 
six- vertex model with reflecting end. To calculate the scalar products and correlation 
functions for the open XXZ model is interesting. 
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